Effects of Chemical Potential on Hadron Masses at Finite Temperature * 
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We study the effects of the chemical potential on the p meson mass at finite temperature. Our preliminary 
results show that some effects are seen in the vicinity of the phase transition point. Although the signal is still 
too noisy to obtain conclusive physical results within limited statistics, the mass susceptibility is consistent with 
zero. 



1. Introduction and Motivation 

QCD Sum Rule analysis has suggested that a 
mass-shift occurs in proportion to the baryon den- 
sity of medium. In nuclear matter, a few to 10 % 
reduction of vector meson masses is expected [Q 
and confirmation of this fact is of current experi- 
mental interest. 

Although the importance of a quantitative in- 
vestigation on this subject is obvious, it seems 
difficult to tackle it by present lattice technol- 
ogy. |0j A new attempt to measure the derivative 
of mass M in terms of the chemical potential, 
dM/dfi at fi = has been made in this work. At 
T = 0, since the density depends cubically on fx, 
dM/dfi vanishes. On the other hand , at finite 
temperature, a linear dependence on the chemical 
potential is expected and we will have a chance 
to see the density effect. 

2. Formulation 

We use 2 flavors of the Wilson fermion for the 
fermionic part of the action. The action is given 
by S = Sc + Sf where Sq and Sf stand for the 
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gauge and the fermionic parts respectively, and 
Sf has the form, 

S F = i'\n)D{n ) m)iP\m). (1) 

i—u : d 

For Sq we use the standard Wilson action. The 
Wilson matrix D(n, m) with chemical potential is 
written asf|,[i 

D(n, m) = 8 n ,rn (2) 
3 
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The hadron correlation function C(t) is given 
by 

C(t) = J dUdipd4>H{t)H\0)e- s /Z (3) 
= 2^coBh(Afi(t-JV t /2)). (4) 

i 

where A\ may be written as 
Ai =A i /{2smh{M i N t /2)). 
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To derive the expression of the mass suscepti- 
bility on the chemical potential, we take a deriva- 
tive of C{t) with respect to the chemical potential 
fi. ^From eq.(0) one obtains 



dC(t) 
dp 



< H(t)H\0) >< ^ > 



(5) 



On the other hand from eq 
r dA t 
dp 

drrij 



dp 

one obtains 



cosh(Mj(t- Nt/2) 



(6) 



dp 



A t {t - Nt/2) sinh(M l (t - Nt/2))}. 



Eqs.fsJ) and @ should be equal to each other. 
Eq.(5) is evaluated by Monte Carlo simulations. 
Eq.(B) contains 4 unknown fitting parameters for 
a given i. A4 and Mi are determined from a fit to 
the hadronic correlation function of eq. (^) . Then 
these results are used for eq.(||) as fixed param- 
eters and, thus the chemical potential dependent 
terms dAi/dp and dMi/dp are extracted from a 
two parameter fit to eq.(^). 

3. Calculation of < H(t)W(0)^- > 

In the present study we aim at obtaining 
the mass susceptibility at p — 0. At p, = 
0, < ^ >= 0. Thus the calculation of 
eq.(pl) reduces to calculate the remaining term 



<if(f)jrt(0)8SE>. 

Now we specify the hadronic operator. For 
the present study we take the p meson operator, 
which is given by, 



H(t) = Y,u(x,t)j v d(x,t). 



(7) 



where 

B, = (1 - l4 )U 4 ( m )e^S^ n+i , (9) 

B 2 = (l + l,)U\{n - 4)e-"* TO>B _ 4 . (10) 

Let us consider 
< H{i)H\Q)K^u{n)B lU {m) > (11) 

where u(m) = tp u (m). This term contains a con- 
nected and disconnected diagrams. We neglect 
the disconnected part, which is unimportant for 
the present study. Eq.(jll]) is written as 

-^7^ _1 (0,0:z)(l~74)C/4(z)x (12) 

x,z 

D- 1 (z + l: x,t) lv D' 1 (x,t : 0,0). 

where the chemical potential p is set to zero. 

Similar calculations are applied for the remain- 
ing terms ( for B2 and for i = d ). 
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We take v — 2 here. The derivative of Sf with 
respect to p is 

dS F 



dp 



(8) 
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Figure 1. Polyakov loop and its susceptibility as 
a function of k. 



4. Simulation 

Simulations are done on a lattice size of 16 x 
8x8x4 at = 5.30. 
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To find a finite temperature transition point, 
the Polyakov loop and its susceptibility are cal- 
culated, as shown in Fig.l The finite temperature 
transition happens around k — 0.158. The chem- 
ical potential effect on hadron masses at /i = 
is expected to start in the vicinity of the tran- 
sition. We take n = 0.153 and k = 0.158 for 
the present study. We use 30 configurations gen- 
erated by the hybrid Monte Carlo algorithm for 
each k. The configurations are separated by 40-50 
trajectories. 

5. Preliminary results 

Figs. 2 and 3 show the derivative of the hadronic 
correlation function C(x) with respect to /x as a 
function of the largest lattice direction: x. The 
hadronic correlation function is measured in the 
x-direction. 
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Figure 2. The derivative of the hadronic corre- 
lator C(x) with respect to // at k = 0.153 as a 
function of x. 



Although the signals within the current statis- 
tics are still too noisy, we fit the data using eqs.Q 
and (||). The mass susceptibility, dM/dfi, is con- 
sistent with zero, while dA/dfi could be finite and 
positive. 
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Figure 3. Same as in Fig. 2 but at k = 0.158. 



6. Discussions 

The effect of chemical potential on the mass 
susceptibility is suppressed by exp(— Mb/T) 
where Mb is the baryon mass. Our results may 
show that Mb at our simulation parameter is too 
heavy to obtain a reasonable signal within lim- 
ited statistics. Our hope is that signal may be 
enhanced if a small quark mass is used. 
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